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Soil  blowing  is  a  serious  problem  in  arid  and  semiarid 
regions,  including  sandy  areas  in  the  United  States.  In  30 
years  of  research,  W.  S.  Chepil  (_j_,  2)  and  his  associates 
(4.  5)  delineated  major  factors  influencing  the  amount 
of  soil  that  wind  will  erode  from  a  given  agricultural 
field.3  They  expressed  these  factors  in  a  wind-erosion 
equation: 

E=f(l',K',C\L',V) 

where  E  is  the  potential  average  annual  soil  loss  for  the 
specified  conditions  of  credibility.  / ;  roughness,  K  ; 
climatic  factor,  C ;  field  length.  L  ;  and  equivalent 
vegetative  cover,  V.  To  facilitate  use  of  this  equation  to 
predict  soil  loss  more  effectively,  and  hence  to  control 
erosion  with  appropriate  practices,  we  wrote  a  Fortran 
IV4  program  WEROS  (Wind  Erosion)  capable  of  solving 
the  equation  quickly  and  accurately.  Besides  predicting 
amount  of  soil  loss  by  wind,  WEROS  (through  high- 
speed calculation)  has  made  it  easier  to  visualize  the 
problem  in  reverse;  that  is,  one  can  specify  (in  the 
program)  the  tolerable  limits  of  soil  erosion,  and  the 
program  then  can  determine  conditions  for  controlling 
soil  loss  within  those  limits.  Thus,  using  WEROS  one  can 
easily  study  the  various  factors  altering  soil  loss. 
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Underlined  numbers  in  parentheses  refer  to  items  listed  in 
List  of  References,  p.  u  . 

Fortran  IV  is  one  of  the  most  widely  used  programing 
languages  for  communicating  with  a  computer. 


This  report  has  two  sections.  The  first  presents  two 
examples,  one  on  predicting  soil  loss  and  the  other  on 
conditions  for  controlling  soil  loss.  Each  example  in- 
cludes information  needed  to  use  WEROS  on  any  similar 
problem.  The  second  section  describes  WEROS  as  a 
computer  program,  along  with  some  relevant  program 
procedures.  To  understand  the  examples,  the  reader 
should  be  familiar  with  information  in  references  (4)  and 
(5). 

PARTI:   HOW  TO  USE  WEROS 

Example  1:  Obtaining  the  soil-erosion  amount  from 
specified  conditions. 

Suppose  we  want  to  know  the  potential  soil  loss 
based  on  existing  conditions.  We  begin  with  the  form  in 
figure  1.  Some  items  in  this  form  are  not  pertinent  to 
this  particular  problem.  For  example,  the  unite  are 
English;  therefore,  the  box  for  metric  units  is  not 
checked.  Similarly,  item  2,  windward  knoll  slope,  does 
not  apply  and  is  left  blank.  In  general,  any  item  on  the 
form  having  an  indicated  default  option  may  be  left 
blank,  resulting  in  use  of  the  default  option.  However, 
all  items  without  options  (Nos.  1,  4,  5,  6,  7)  must 
contain  some  information  for  WEROS  to  obtain  an 
answer.  Information  that  does  not  apply  can  be  left 
blank,  implying  a  value  of  zero  for  numerical  items,  and 
"does-not-apply"  for  checked  items.  Item  8  (left  blank 
here)  pertains  only  when  form  2  (fig.  2)  follows.  Form  2 
specifies  the  information  that  WEROS  uses  to  search  for 
some  condition  to  obtain  the  prespecified  potential  soil 
loss;  it  will  be  discussed  in  example  2. 

Data  on  form  1  are  punched  on  four  cards  in  fields, 
each  10  columns  wide.  Thus,  each  card  will  contain  eight 
items  located  anywhere  in  the  field.  Each  item  will  be 
punched  with  an  explicit  scale  or  decimal  point  or  left 
blank  as  applicable.  Figure  3  shows  the  data  cards 
prepared  to  represent  the  specified  data  for  this  prob- 
lem. The  numerical  values  are  represented  explicitly;  the 


FORM  1  FOR  COMPUTER  SOLUTION  OF  WIND  EROSION  EQUATION 

If  Che  Information  for  this  form  has  metric  units,  check  /~7. 

1.  Percentage  of  soil  aggregates  greater  than  0.84  mm.  /O  % . 

2.  Windward  knoll  slope,  less  500  feet Z. 

a.  From  top  of  knoll  f~J 

b.  From  upper  1/3   /~7 
Default  Is  -  1 

3.  Fill  in  the  appropriate  (not  more  than  one  entry): 

a.  Ridge  height (in.  or  cm.) 

Ridge  spacing (in.  or  cm.) 

b.  Ridge  roughness  factor  (dimensionless) 

c.  Field  is  smooth  ZZ7  or  rough  /"»7l   Default  K'  -  1. 

4.  The  climatic  factor  is  /  QQ    Z. 

5.  Fill  in  the  following: 

a.  The  angle  of  deviation  of  prevailing  wind  erosion  direction  from 
right  angles  to  field  strip  is  ft        degrees. 

b.  Preponderance  of  wind  erosion  forces  in  prevailing  wind  erosion 
direction  is  '2 

c.  Height  of  barrier  (ft.  or  m.) 

d.  Field  width  /yp  (ft.  or  m.).   Default  -  5,000  feet. 

6.  Quantity  of  vegetative  cover  jfoQO  (lbs. /acre  or  kg. /ha.) 

7.  Type  of  vegetative  cover  (check  all  correct  responses): 

a.  Small  grain  stubble  /~~7 

1.  Flat   £7 

2.  Standing  O 

b.  Small  grain  crops  in  seedling  or  s tooling  stage  /~~7 

1 .  Furrow  /~7 

2.  Smooth   O 

c.  Sorghum  stubble  /  *•/ 

1.  Flat  d 

2.  Standing  QET 

3.  Height  IZ.     (in-  or  cm.) 

8.  Check  if  Form  2  follows   £7. 

Figure  1.— Form  1  containing  the  specified  information  for  example  1. 


FORM  2  FOR  COMPUTER  SOLUTION  OF  WIND  EROSION  EQUATION 


If  the  information  for  this  form  has  metric  units,  check  CU 

1.   What  is  tolerable  amount  of  wind  erosion  (tons/acre/year  or  metric 
tons /ha. /year) ? 


Lower  limit   >"""    Upper  limit  £T        Increment 


2.  Check  the  blank,  opposite  the  variable  under  investigation.  Check 
the  box  opposite  the  variable  specified  and  specify  this  variable 
by  entering  appropriate  values. 

'    [~1     a.   Soil  cloddiness  (percentage  of  aggregates  greater  than 
0.84  mm.) 

Lower  limit Upper  limit Increment 

l~l     b.   Soil  ridge  roughness  factor  (dimensionless) 

Lower  limit Upper  limit Increment 


/Z7  c.   Quantity  of  vegetative  cover  (lbs. /acre  or  kg. /ha.) 

Lower  limit Upper  limit Increment 

_  EJ     d.   Barrier  height  (ft.  or  m.) 

Lower  limit Upper  limit Increment 

_  O  e.   Field  width  (ft.  of  m.) 

Lower  limit Upper  limit Increment 

Hn     f .   Climatic  factor  (percent) 

Lower  limit  /Q      Upper  limit  /QO    Increment  /O 

The  program  searches  the  conditions  necessary  to  limit  potential  soil 
loss  to  amount  specified  in  item  1,  +  %.   Default  -  5Z. 

Figure  2.— Form  containing  the  specified  information  for  example  2. 


^~ 


A 


.5 
« 

c 
o 

u 

e 

o 

'43 


g) 

■43 

c 


■3 


CO 


CN 


.V 


checks  on  the  form  are  represented  on  the  cards  by  a 
"1."  Item  number  or  numbers  given  below  each 
10-column  field  correspond  to  the  applicable  numbers 
on  form  1 . 

Data  cards  are  then  placed  at  the  back  of  the 
program,  and  with  the  proper  control  cards  (those 
required  by  the  particular  computer  installation)  in 
place,  the  program  can  be  executed. 

The  output  from  W  EROS,  shown  in  figure  4,  consists 
of  three  basic  items:  (1)  A  brief  copy  of  the  input  data; 
(2)  the  calculated  answer;  and  (3)  identification  informa- 
tion (that  is,  example  1).  The  input  data  as  printed  out 
agrees  with  the  conditions  we  specified  initially  for  this 
problem,  and  the  answer  is  given  as  28  tons/ acre/year  or 
63  metric  tons/hectare/year. 

To  understand  the  mechanics  of  this  procedure;  that 
is,  use  of  data  cards  and  interpreting  information  as  it  is 
reproduced  by  the  computer,  the  following  information 
is  pertinent.  As  shown  in  figure  4,  under  "INPUT,"  the 
units,  when  presented,  are  given  in  the  English  system. 
However,  in  no  way  is  the  card  input  limited  to  English 
units;  in  fact,  if  the  box  denoting  metric  units  is 
checked,  metric  units  must  be  given  for  all  items  or 
answers  will  be  incorrect.  Also  the  number  0.  appearing 
in  figure  3  beside  various  entries  denotes  either  that  the 
item  is  not  applicable,  or  that  its  value  is  0..  The  No.  1. 
has  the  same  function;  that  is,  it  means  that  the  item 
marked  is  applicable,  or  that  its  value  is  1..  Consider,  for 
example,  the  entry  "CLIMATIC  FACTOR."  The  value 


given  is  1.  (fig.  4),  which  means  100  percent,  and  in  this 
case,  it  is  the  literal  value  of  the  climatic  factor  used. 
However,  under  ridge  roughness  (item  3,  fig.  1)  a  1. 
appears  next  to  the  entry  "OR  ROUGH,"  which  signifies 
that  the  field  under  consideration  is  rough.  Another 
anomaly  appears  as  part  of  the  last  entry.  Under 
"SORGHUM  STUBBLE"  and  next  to  "OR  STAND- 
ING" is  the  number  pair,  "1.  12.".  The  1.  again  tells  us 
that  we  are  concerned  with  standing  sorghum  stubble; 
12.  denotes,  in  inches,  the  height  of  the  stubble. 

Returning  to  form  1  (fig.  1)  and  fig.  3,  all  percentage 
data  are  to  be  punched  on  the  cards  as  percentages,  not 
as  decimals.  For  example,  in  figure  3  the  first  card 
contains  an  entry  for  the  percentage  of  soil  aggregates 
larger  than  0.84  mm.;  the  value  given,  10  percent,  is 
punched  10.  on  the  data  card.  Those  entries  not 
applying  or  having  a  value  of  0  can  be  represented  by  a 
blank  or  empty  field  on  the  data  card.  Last,  nothing 
other  than  numbers  and  their  accompanying  decimal 
points  should  appear  on  the  data  cards. 

Example  2:  Obtaining  a  specified  condition  to  con- 
trol erosion  when  tolerable  soil  loss  has  been  specified. 

Consider  a  problem  with  these  specifications: 

1 .  From  a  specified  field,  we  wish  to  allow  a  soil  loss 
of  approximately  5  tons/acre/year. 

2.  The  field  is  smooth. 

3.  The  preponderance  of  wind-erosion  forces  in  the 
prevailing  wind-erosion  direction  is  2,  while  the  angle  of 
the  prevailing  wind-erosion  direction  from  right  angles  to 
the  field  strip  is  20  degrees. 


OUTPUT 


THE  ANSWER  IS 


28.1SS9 


TONS/ACRE/YEAR 


OR 


63.2129         MET    TONS/HA/YEAR 


INPUT 

PERCENTAGE  CF  SOIL  AGGREGATES  10. 0 

WINDWARD  KNCLL  SLOPE   0.0 

A.    0.  R.    0. 

RI0GE  HT  0.C0   IN. 

RIDGE  ROUGHNESS   0.0 

SMOOTH    0.       OR  ROUGH    1. 

CLIMATIC  FACTOR    1.0 

WIND  ANGLE   0.0 

WIND  FORCE    2.00 

HFIGH  OF  eARRIER  (FT)      0.0 

FIELD  WIDTH  (FTI   1320.0 

VEGETATIVE  COVER  (LBS/ACRE)   1600.0 

SMALL  GRAIN  STURBL^    0. 

FLAT    0.      OR  STANDING    0. 
SMALL  GRAIN  SEEDLING  OR  STOOLING    0. 

FURROW    0.      OR  SMOOTH    0. 
SORGHUM  STueBLE 

FLAT    0.      OR  STANDING    1.   12. 


KIDGE  SPACING    O.JO   IN. 


CASE  NUMBER  = 


Figure  4.— Output  of  WEROS  after  processing  the  information  of  example  1,  as  specified  on  the 

cards  shown  in  figure  3. 


4.  The  field  is  1,480  feet  wide. 

5.  It  has  no  vegetative  cover. 

Using  those  specifications,  suppose  we  investigate  the 
effect  of  climatic  factor  on  the  percentage  of  soil 
aggregates  greater  than  0.84  mm.  required  to  limit  soil 
loss  to  a  specified  amount.  That  is,  for  different  climatic 
factors  what  percentage  of  soil  aggregates  greater  than 
0.84  mm.  is  needed  to  keep  wind  erosion  to  the  level  of 
5  tons/acre/year?  Let  us  examine  climatic  factors,  with 
an  increment  value  of  10,  ranging  from  a  low  of  10  to  a 
high  of  100. 

Form  1  can  be  completed  as  shown  in  figure  5; 
however,  this  time  item  8  must  he  checked.  That 
signifies  to  WEROS  that  what  follows  is  a  data  card 
applying  to  the  search  part  of  the  program  and  supplying 
the  necessary  information  to  ohtain  the  required  answer. 

With  item  8  of  form  1  checked,  we  now  complete 
form  2  as  shown  in  figure  2.  Again  the  units  are  in 
English,  so  the  first  hox  on  the  form  is  left  blank.  As 
specified,  we  are  interested  in  a  soil  loss  of  approxi- 
mately 5  tons;  therefore,  in  item  1,  we  will  supply  5.  for 
both  the  lower  and  the  upper  limits,  leaving  the 
increment  blank. 

In  item  2  we  are  interested  in  the  effect  of  climatic 
factor  on  soil  cloddiness,  so  we  check  the  blank  in  part  a 
and  the  box  in  part  f,  here  supplying  the  lower  (10.)  and 
upper  (100.)  limits  of  the  climatic  factor  and  its 
increment  (10.). 

Item  3  allows  us  to  adjust  the  tolerance  on  soil-loss 
limits  recorded  in  item  2.  That  is,  if  item  3  is  left  blank, 
the  program  will  accept  an  answer  for  soil  cloddiness 
that  under  the  specified  climatic  condition  will  be  in  the 
range  of  5  tons/acre/year  ±  5  percent,  or  4.75  to  5.25 
tons/acre/year,  but  if  one  wishes  an  answer  to  any  other 
degree  of  accuracy  that  value  can  be  placed  in  item  3. 
For  this  problem,  the  5  percent  limit  is  deemed 
acceptable. 

Figure  6  shows  the  data  cards  prepared  for  this 
problem.  The  first  four  cards,  representing  information 
on  form  1,  have  been  punched  as  specified  in  the 
previous  example.  The  fifth  card,  representing  informa- 
tion on  form  2,  contains  at  most  10  punched  numbers, 
each  number  being  allowed  five  columns.  The  first  five 
columns  contain  either  a  blank  or  the  right  justified 
number  "1,"  depending  on  whether  English  or  metric 
units  are  used.  The  next  three  sets  of  five  columns 
(columns  6-20)  contain  the  lower,  upper,  and  incre- 
mental values  for  soil  loss,  respectively;  each  of  the  three 
numbers  includes  a  decimal  point  and  can  be  located 
anywhere  in  its  designated  5-column  field.  The  next  two 
fields  (columns  21-30)  contain  information  for  item  2; 
the  first  number  is,  for  this  example,  a  "1,"  right 
justified,  denoting  that  the  first  entry  of  item  2  is  the 


variable  we  wish  to  determine.  The  second  field  of  5 
columns  contains  a  "6,"  right  justified,  signifying  that 
entry  "f"  (sixth  entry,  climatic  factor)  is  the  one  that 
will  vary.  The  next  three  fields  (columns  31-45)  contain 
the  numbers  representing  lower,  upper,  and  incremental 
values  of  the  climatic  factor  as  specified;  the  values 
should  include  a  decimal  point  and  may  be  located 
anywhere  in  their  respective  fields.  If  no  decimal  is 
given,  the  program  automatically  places  a  decimal  one 
digit  to  the  left  of  the  last  digit  in  the  field. 

To  instruct  the  computer  to  print  out  each  iteration 
of  the  search,  add  1,000  to  the  increment  value.  The 
increment  value  (the  "10"  of  fig.  2),  then,  is  punched  as 
10,100  as  shown  in  figure  6.  The  program  interprets  this 
as  1,010.0  and  will  extract  the  increment  value.  It  sets  a 
switch  for  this  case  and  prints  out  each  iteration  of  the 
search  as  shown  in  figure  7. 

The  output  from  this  problem  (shown  in  fig.  7) 
consists  of  the  important  input  parameters  and  their 
values  as  described  in  the  last  section.  The  identification 
"CASE  NUMBER  =  2"  is  included  to  indicate  this  was 
the  second  problem  run.  The  next  row  shown  pertains  to 
data  on  form  2;  that  is,  we  are  looking  for  1  (soil 
cloddiness)  using  6  (climatic  factor)  where  the  values  of 
6  range  from  a  low  of  0.1  to  a  high  of  1.  by  increments 
of  0.1  or  in  percentages  10,  100,  10,  respectively.  The 
next  item  indicates  amount  of  soil  loss  specified  for  the 
conditions  specified  in  the  preceding  line— in  this  case,  5 
tons/acre/year. 

Because  the  increment  value  of  0.1  was  punched  on 
the  card  as  10,001,  the  next  five  rows  are  printed  out  to 
show  the  search  number  (that  is,  each  try  of  the  program 
to  find  the  proper  soil-cloddiness  value).  The  second 
entry  on  each  of  these  five  lines  contains  information 
that  prevents  the  program  from  continuing  indefinitely 
without  converging  to  an  answer.  In  general,  if  the  next 
guess  does  not  improve  the  answer,  then,  the  "improve- 
ment-counter" is  incremented.  If  five  consecutive  non- 
improvements  are  recorded,  the  search  will  terminate. 
The  program  will  then  proceed  to  the  next  case,  be  it  a 
new  set  of  data  or  just  the  increment  of  a  current 
parameter.  This  counter  is  used  primarily  when  a  set  of 
conditions  is  given  in  which  no  answer  is  possible,  and 
with  our  example,  that  could  occur  if  the  percentage  of 
soil  cloddiness  needed  to  control  soil  erosion  happened 
to  be  outside  the  range  of  the  soil-cloddiness  table. 

The  search  bounds,  given  next  in  these  five  rows, 
specify  the  lower  and  upper  limits  of  the  interval  in 
which  the  value  of  cloddiness  is  being  sought.  The  last 
two  items  of  each  row  specify  the  result  of  soil  loss  for 
the  value  given.  Thus,  we  read  the  first  of  these  five  rows 
of  search  information  in  figure  7  as:  It  is  search  No.  1  in 
the  interval  2.  to  80.  (soil-cloddiness  range),  and  WEROS 


FORM  1  FOR  COMPUTER  SOLUTION  OF  WIND-EROSION  EQUATION 

If  the  information  for  this  form  has  metric  units,  check  /~7. 

1.  Percentage  of  soil  aggregates  greater  than  0.84  mm.  . 

2.  Windward  knoll  slope,  less  500  feet %. 

a.  From  top  of  knoll  f~l 

b.  From  upper  1/3    f~l 
Default  I8  -  1 

3.  Fill  in  the  appropriate  (not  more  than  one  entry): 

a.  Ridge  height (in.  or  cm.) 

Ridge  spacing  (in.  or  cm.) 

b.  Ridge  roughness  factor  (dimensionless) 

c.  Field  is  smooth  HT  or  rough  CJ •   Default  K'  -  1. 

4.  The  climatic  factor  is  %. 

5.  Fill  in  the  following: 

a.  The  angle  of  deviation  of  prevailing  wind  erosion  direction  from 
right  angles  to  field  strip  is  "2Q        degrees. 

b.  Preponderance  of  wind  erosion  forces  in  prevailing  wind  erosion 
direction  is  ? 

c.  Height  of  barrier (ft.  or  m.) 

d.  Field  width  Jf/SO   (ft.  or  m.).   Default  -  5,000  feet. 

6.  Quantity  of  vegetative  cover (lbs/acre  or  kg. /ha.) 

7.  Type  of  vegetative  cover  (check  all  correct  *•  .ponses): 

a.  Small  grain  stubble  CJf 

1.  Flat  O 

2.  Standing  O 

b.  Small  grain  crops  in  seedling  or  s tooling  stage   /~~7 

1 .  Furrow  /~7 

2.  Smooth  O 

c.  Sorghum  stubble   /~~7 

1.  Flat  O 

2.  Standing  EJ 

3.  Height  (in.  or  cm.) 

8.  Check  if  Form  2  follows  CB^. 

Figure  5.— Form  1  containing  the  specified  information  for  example  2. 
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INPUT 

PERCENTAGE  CF  SOIL  AGGREGATES   0.0 

WINOWARO  KNCLL  SLOPE   0.0 

A.    0.  P.    0. 

RIOGE  HT  O.CO   IN. 

RIDGE  ROUGHNESS   0.0 

SMOOTH    1.      OR  ROUGH 

CLIMATIC  FACTOR    0.0 

WIND  ANGLE  20.0 

WIND  FORCE    2.00 

HEIGH  OF  BARRIER  ( FT  I 

FIELD  WIDTH  (FT)   1480.0 

VEGETATIVE  COVER  (LBS/ACRE) 

SMALL  GRAIN  STUBBLE    0. 

FLAT    0.      OR  STANDING    0. 
SMALL  GRAIN  SEEDLING  OR  STCOLING    0. 

FURROW    0.      OR  SMOOTH    0. 
SORGHUM  STUBBLE 

FLAT    0.      OR  STANDING    0.    0 


RIDGE  SPACING    0.00   IN. 


0. 


0.0 


0.0 


CASE  NUMBER  = 


LCOKING  FOR 
NEXT  CASE  FOR  E= 


SEARCH  NO. 
SEARCH  NO. 
SEARCH  NO. 
SEARCH  NO. 


1  NO. 

2  NO. 

3  NO. 

4  NO. 


SEARCH  NO.     5   NO. 
NUMBER  OF  SEARCHES 
GIVEN  NO     6 
SOUGHT   VALUE 


SEARCH  NO.  1 
SEARCH  NO.  2 
SEARCH  NO.     3 


SEARCH  NO.     4 
NUMBER  CF  SEARCHES 
GIVEN  NC     6 
SOUGHT   VALUE     1 


1 

NO. 
NO. 
NO. 

NO. 


1   US1 
5.0 

FROM 

FROM 

FROM 

FROM 

FROM 
5 

FROM 
FROM 
FROM 

FROM 

4 


NG 

LAST 
LAST 
LAST 
LAST 
LAST 

LAST 
LAST 
LAST 
LAST 


6   VALUES  OF: 


LOWER 


0.10   UPPER 


I. 00   AT  INCREMENTS  OF; 


o.l. 


IMPROV. 
IMPROV. 
IMPROV. 
IMPROV. 
IMPROV. 


0  SEARCH 

0  StARCH 

0  SEARCH 

0  SEARCH 

0  SEARCH 


BOUNDS 
BOUNDS 
BOUNDS 
BOUNDS 
BOUNDS 


2.0000  80 

2.0C00  8C 

2.0000  41 

21.5000  41 

31.2500  41 


♦♦RESULT** 
IMPROV.     0   SEARCH 


IMPROV. 
IMPROV. 
IMPROV. 


0  SEARCH 
0  SEARCH 
0   SEARCH 


SEARCH  NO. 


NO.  FROM  LAST 


♦♦RESULT** 
IMPROV.      0   SEARCH 


36.13  (ENGLISH  UNITS) 
BOUNDS  2.0000  80. 
BOUNDS  2.0000  80, 
BOUNDS  41.0000  80, 
BOUNDS    41.0000    60, 

50.75  (ENGLISH  UNITS) 
BOUNDS     2.0000    80 


0000   ANS  24,9999   FOX    2.000 

0000   ANS  4.2425   FOK   41.000 

,0000   ANS  8.8926   FOR   21.5u0 

,0000   ANS  5.4632   FOk   31.25u 

0000   ANS  4,7750   FOK   36.125 

36.  13  (METRIC  UNITS) 

,0000   ANS  49.9288   FOk    2.000 

,0000   ANS  8.5636   FOR   41.000 

,0000   ANS  2.5906   FUR   60.500 

,5000   ANS  4,d569   FOR   50.750 

50.75  (KtTKlC  UNITS) 

,0000   ANS  74.9020   FOR    2.000 


Figure  7.— Output  of  WEROS  for  example  2. 


has  tried  a  value  of  2.,  yielding  a  soil  loss  of  24.999  or 
25.  tons/acre/year.  That  is  far  in  excess  of  the  5 
tons/acre/year  specified.  Therefore,  the  process  is  con- 
tinued until  (at  the  completion  of  the  fifth  search) 
WEROS  has  obtained  a  soil-cloddiness  value  of  36.125 
percent,  which  yields  a  soil  loss  of  4.775;  that  is,  within 
the  5  percent  range  of  5  tons/acre/year  for  a  climatic 
factor  of  10. 

The  next  three  lines  specify  the  given  and  sought 
value  as  well  as  the  solution.  In  this  case,  the  answer  for 
English  and  metric  units  is  identical  because  the  answer 
is  a  percentage  without  units.  However,  when  the  answer 


has  units,  then  both  the  value  in  metric  units  and  the 
value  in  English  units  will  be  given. 

WEROS  then  continues  this  case  by  incrementing  the 
climatic  factor  of  10  by  10  to  get  20  and  proceeds  to 
compute  a  new  soil-cloddiness  lac  tor,  which  it  does  in 
four  searches.  That  value  is  given  as  50.75  percent. 

Item  1  on  form  2  can  specify  three  things.  One 
situation  investigated  was  the  effect  of  varying  climatic 
factor  on  soil  cloddiness  with  soil  loss  specified  at  5 
tons/acre/year.  Suppose  now,  we  investigate  a  pair  of 
factors  at  some  specified  soil  los;-:  for  example,  the 
effect  of  climatic  factor  and  soil  cloddiness.  We  can  do 


that  bv  specifying  llie  limits  and  an  increment  value.  For 
example,  suppose  we  apply  the  last  problem  to  values  of 
soil  loss  from  5  lo  25  tons/acre/year  at  increments  of  5 
tons.  Tliis  information  should  be  included  on  the  input 
card,  and  the  output  should  contain  (for  the  above 
problem)  10  answers  for  3  Ions.  10  for  10  tons/acre/ 
year,  and  >o  forth  lor  a  total  of  50  answers. 

Then  we  interpret  item  I  of  form  2.  Suppose  we  are 
now  interested  in  the  effect  of  climatic  factor  on  soil 
cloddiness  (again  we  could  have  selected  any  other  pair) 
when  the  soil  loss  is  specified  onl)  h\  a  low  and  a  high 
value.  I' or  example,  we  could  allow  soil  loss  lo  be 
specified  in  the  range  5  lo  10  tons/ acre/ \  ear;  then  5 
would  he  included  in  the  "lower  and  10  in  the  "upper" 
and  the  increment  left  blank.  Thus,  the  example 
explained  above  is  just  a  special  case  of  this  more  general 
abililj . 

The  last  items  of  importance  to  the  users  of  WEROS 
arc:  (I)  Form  I  and  form  2  need  not  specif)  the  same 
type  of  units,  and  (2)  when  both  forms  are  used,  any 
information  contained  on  the  cards  representing  form  1, 
which  pertains  to  or  duplicates  information  given  on  the 
card  representing  form  2.  will  be  disregarded.  That  is.  if 
in  example  2,  the  cards  from  form  I  bad  contained  a 
value  of  2.  for  the  climatic  factor,  this  value  would  have 
been  ignored  for  the  computations  specified  by  the  card 
representing  form  2,  which  contains  conflicting  specifi- 
cations lor  the  climatic  factor. 


PART  2:   THE  PROGRAM  WEROS 

WEROS  is  written  in  a  generally  accepted  implemen- 
tation—an  independent  version  of  FORTRAN  IV.  For 
the  most  part.  WEROS  follows  the  procedure  outlined 
by  Skidmorc  and  others  (.})  to  obtain  a  solution  when 
the  problem  is  similar  to  that  of  example  I  of  the 
previous  part.  However,  the  solution  procedure  is  some- 
what more  complex  when  the  problem  is  similar  to  that 
given  in  example  2  of  Part  I . 

WEROS  consists  of  1  I  subroutines,  each  performing  a 
major  function  in  the  solution  procedure.  Thus,  the 
program  can  be  modified  and  tailored  easily  lo  fit 
individual  requirements.  This  modularity  also  contrib- 
utes lo  understanding  of  functions  performed  in  each 
routine.  Figure  8  shows  the  sequence  of  routines  as  the) 
are  used  to  solve  problems  similar  to  examples  1  and  2 
of  the  previous  part,  respectively. 

The  only  significant  difference  between  the  two 
sequences  (fig.  8)  is  the  use  of  MAIN  1  when  a  problem 
requiring  a  search  is  specified. 

A  brief  explanation  of  each  subroutine  follows. 


A.  MAIN 

MAIN  is  the  controlling  routine,  and  as  shown  in 
figure  8,  it  will  call  for  the  routine  IPDT;  it  may  call  for 
the  routine  COMPL'T  or  MAIN1.  depending  on  the 
specification  of  the  problem.  IPUT  is  called  to  obtain 
data  from  the  necessary  tables  and  figures  as  contained 
in  (4).  Data  representing  information  on  form  1  are  read 
in  next,  and  if  these  data  are  in  metric  units,  they  will  be 
converted  to  the  English  equivalents.  Thus,  all  computa- 
tions in  WEROS  are  done  with  Fnglish  units,  which 
simplifies  programing. 

II  method  I  (that  is  the  procedure  required  to  solve 
problems  similar  to  example  I)  is  required,  MAIN  calls 
COMPl  T.  which  will  obtain  the  answer.  The  answer  is 
[irinted  out:  then,  the  input  data  are  printed,  completing 
the  required  task.  If.  however,  method  2  (lhal  is.  the 
procedure  required  to  soke  problems  similar  lo  example 
2)  i>  specified,  then  MAIN  prints  out  the  input  data  and 
calls  MAIN  I  to  complete  the  problem.  As  each  specified 
problem  is  completed,  another  set  of  data  cards, 
representing  form  I,  is  read  in.  If  no  such  cards  arc 
found,  WEROS  is  terminated  with  an  error  message  due 
to  an  end-of-file  condition.  This  termination  procedure 
could  easily  be  altered  lo  meet  the  specifications  of  an 
individual  system  by  inserting  statements  lo  check  for  a 
special  terminal  card. 

B.  INPUT 

The  INPUT  routine  consists  of  read  statements, 
which  read  from  cards,  the  data  representing  the  tables 
and  graphs  used  in  calculating  an  answer.  The  data  arc 
stored  in  seven  arrays  as  discrete  values  from  the  eurws 
they  represent.  Oilier  values  that  may  be  required  are 
interpolated  from  surrounding  values.  We  used  values 
obtained  from  graphs  based  on  the  criteria  that  in  the 
areas  where  the  slopes  of  the  curves  exhibited  the  most 
rapid  change,  the  values  were  taken  closer  together, 
while  from  the  remaining  areas  of  the  graphs  fewer 
values  were  used.  Selecting  the  points  in  this  manner 
permitted  use  of  linear  interpolation  with  a  relatively 
small,  but  predictable  error. 

C.  COMPUT 

This  routine  controls  the  remaining  routines  and  the 
sequence  in  which  they  are  executed  to  obtain  an  answer 
for  the  specified  data. 

D.  IPRIME 

This  routine  calculates  the  value  of  /,  which  is  soil 
erodibility,   from   the  input  data  for  the  percentage  of 
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A. 


'NEXT  CASE* 


COMPI  T •-  FIN 

COMPUT  > 
MAIN  ^— *► 'END  OF  DATA' 


'INTER 
'INTER 
[INTER 

INTER 


'NEXT  CASE' 


MAIN 


MAINF 
+-  'END  OF  DATA' 


7  INTER 

Winter 

^  INTER 
INTER 


Figure  8.— Sequence  of  subroutines  to  calculate  an  answer  to  problems  similar  to  examples  A  and 
B,  respectively.  The  broken  lines  represent  possible  use;  the  solid  Unes.  definite  use. 


soil  fractions  larger  than  0.84  mm.  Then,  using  tin 
information  on  the  percentage  of  windward  knoll  slope, 
the  value  of  1  (soil-erodibility  index)  is  calculated.  If  the 
percentage  of  windward  knoll  slope  is  given,  then 
1PR1ME  will  use  INTER  to  obtain  the  potential  soil  loss 
for  that  value.  If  no  percentage  of  windward  knoll  slope 
is  specified,  the  value  of  1  will  be  used. 


E.  KPRIME 

This  routine  determines  the  soil-ridge  roughness 
factor  K  from  data  specified  for  this  item  on  form  1 .  If 
a  measurement   for  ridge  spacing  and   height  is  given. 


KPRIME  will  use  INTER  to  interpolate  a  value  for  the 
soil-ridge  roughness  factor  K .  If  the  field  is  simply 
smooth  or  rough,  then  KPRIME  will  return  the  value  1. 
or  0.5.  respectively.  If  nothing  is  specified.  KPRIME  will 
assume  a  smooth  field  and  return  the  value  I . 

F.   LPR1ME 

I  hi  routine  calculates  the  median  unsheltered  travel 
distan< -;•  across  field,  L  .  It  involves  items  listed  under  5 
on  form  1.  In  the  worst  ease,  three  interpolations  are 
required  to  obtain  the  median  travel  distance  across 
field,  from  which  LPRJME  subtracts  10  times  an) 
barrier  height  to  obtain  unsheltered  distance. 
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G.  FLRM 

This  routine  uses  the  value  oi  L  ,  E2,  E3  to  calculate 
a  value  of  E4  where 

E2  -  /'  K' 

E3  =/>2  C  (climatic  factor). 

£4  is  the  soil  loss  from  a  field  having  no  vegetative 
cover.  En  is  expressed  as 

E,=l'K'C'f(L')-E3f(L'). 

To  calculate  E4 ,  a  series  of  interpolations,  or  manipu- 
lations, of  two  figures  are  required.  The  manipulation,  as 
specified  in  (4,  5),  is  impossible  to  do  with  a  computer. 
Therefore,  an  equivalent  procedure,  explained  in  (3),  is 
used.  Use  of  the  new  approach  requires  several  interpola- 
tions that  cannot  he  included  conveniently  in  routine 
INTER.  Therefore,  they  are  included  within  FLRM,  and 
they  follow  the  procedure  outlined  in  (3). 

H.  VEG 

This  routine  uses  a  set  of  nine  graphs  to  determine 
the  equivalent  pounds  of  vegetative  cover  on  a  field. 
Here  again,  depending  on  conditions  specified  by  form 
I .  one  or  three  interpolations  could  be  necessary.  Using 
the  specified  data,  VEG  will  use  INTER  to  obtain  the 
value  V. 

1.    M1N 

This  routine  uses  the  values  of  V  (equivalent  pounds 
per  acre  of  vegetative  cover)  and  E4  to  determine  the 
answer  of  soil  loss.  The  family  of  curves  given  in  (4,  fig. 
24)  and  in  (5,  fig.  10)  shows  a  sequence  of  polynomials 
of  order  one  when  they  are  graphed  on  a  log-log  scale. 
Being  straight  lines,  they  can  be  stored  as  a  pair  (x,6), 
where  x  is  logj  0  of  the  abscissa  intersection  value  and  6 
is  the  acute  angle  between  the  abscissa  and  the  line.  FIN 
takes  the  logarithm  of  E4  and  finds  the  difference 
between  this  value  and  the  point  of  intersection  of 
straight  line  with  abscissa.  That  difference,  which  repre- 
sents the  adjacent  side  of  a  right  triangle  whose 
hypotenuse  is  equivalent  vegetative  cover,  and  angle  6 
are  used  to  determine  the  ordinate,  which  is  the  log  of 
sought  answer.  If  the  value  of  V  given  to  FIN  is  not  one 
of  the  exact  values  of  V  on  the  curves,  interpolation  is 
required  and  two  ordinate  values  are  determined  from 
the  formula 

yy  =  xytanQi 
Yi  =  x2tanQ2 


where  0  j  and  62  are  the  angles  of  the  preceding  and 
following  lines  corresponding  to  the  surrounding  value 
of  V,  and  x\  and  x2  are  the  abscissa  values  corre- 
sponding to  6 1  and  62.  Using  the  antilog  of  ji  andy2. 
an  ordinate  value  y  can  be  calculated  by 


y-y\  ~v  _yy'(yi  -y^ 


where  y\   and  v'2  are  the  antilogs  of  yj  and  y2  and  Vx 
and  V2  arc  lines  preceding  and  following  V,  respectively. 

J.  INTER 

INTER  is  a  general  purpose  linear-interpolation  rou- 
tine; it  can  interpolate  on  either  abscissa  or  ordinate 
values,  depending  on  the  value  of  an  input  parameter. 

K.  MAINl 

This  routine  is  used  when  method  2  is  specified.  Like 
MAIN,  it  has  several  functions.  The  first  is  to  read  in  the 
data  card  representing  form  2  and  based  on  what  is  con- 
tained on  this  card,  MA1N1  defines  the  problem  to  be 
solved.  The  process  involves  these  steps: 

1.  Define  the  upper  and  lower  limits  allowed  for  soil 
loss. 

2.  Determine  the  bounds  of  the  interval  of  definition 
for  the  variable  being  investigated  and  assign  thai 
variable  a  value  equal  to  its  lower  limit. 

3.  Assign  the  initial  value  of  the  specified  variable  to 
its  proper  name  with  WEROS. 

4.  With  this  information,  call  COMPUT,  which  pro- 
ceeds with  the  calculation  already  specified. 

5.  After  an  E  value  or  answer  is  calculated,  return  to 
MUNI  and  check  to  see  if  it  is  within  the  necessary 
hounds.  If  it  is,  then  continue  with  the  next  specified 
task.  If  it  is  not,  adjust  the  value  of  the  variable  to  a  new 
value  and  call  COMPUT  to  continue  processing. 

The  method  of  search  across  the  interval  of  definition 
for  the  variable  being  investigated  is  the  well-known 
halving  technique.  That  is,  the  first  two  points  are  the 
interval  bounds;  the  third  is  the  center  point.  The  three 
points  are  used  to  select  a  new  interval,  which  is  one-half 
the  old,  and  a  new  value,  which  is  the  midpoint  of  the 
new  interval.  Thus  in  re  >  3  searches,  WEROS  is  working 
with  an  interval  of  Vz  re-2  times  the  length  of  the  original 
interval.  The  number  of  searches  generally  ranges  from 
four  to  six. 

As  explained  in  example  2  of  Part  1 ,  several  values  of 
E  and  several  values  of  the  specified  variable  could  be 
given  by  stating  a  lower  and  an  upper  limit  with  an 
increment  value.  WEROS  uses  the  values  of  £  as  the 
outer  loop  and  the  values  of  the  specified  variable  as  the 
inner  loop. 
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CONCLUSION 


Although  the  brief  descriptions  included  here  do  not 
give  all  details  of  WEROS.  at  least  function  and  role  of 
each  routine  are  specified  regarding  the  problem-solving 
process  embodied  in  the  program.  If  source  deck  is 
desired,  it  can  be  obtained  by  writing  to  the  authors.  In 


requesting  a  source  deck,  include  the  type  of  machine 
facility  that  is  available  at  the  requesting  facility  to 
insure  as  little  trouble  as  possible  in  getting  WEROS  to 
function  correctly. 
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